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Abstract 

Let B H and B H be two independent, d-dimensional fractional Brow- 
nian motions with Hurst parameter H £ (0, 1) . Assume d > 2. We prove 
that the intersection local time of B H and B H 



exists in L 2 if and only if Hd < 2. 
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1 Introduction 

We consider two independent fractional Brownian motions on M. d , d > 2, with the 
same Hurst parameter H £ (0, 1) . This means that we have two d-dimensional 
independent centered Gaussian processes B H = {B^ , t > 0} and B H = {B^ 1 , t > 
0} with covariance structure given by 




^[Bf^Bf'i 



E[B?< i B?>i]=5 ij R H (M), 



where i, j = 1, d, s, t > and 



R H (s,t) = - (t 



2H 



+ S 



2H 




The object of study in this paper will be the intersection local time of B and 
B H , which is formally defined as 

I(B H ,B H )= / / So{Bf - Bf)dsdt, 
Jo Jo 

where Sq (x) is the Dirac delta function. It is a measure of the amount of time 
that the trajectories of the two processes, B H and B H , intersect on the time 
interval [0, T] . As we pointed out before, this definition is only formal. In order 
to give a rigorous meaning to I(B H , B H ) we approximate the Dirac function by 
the heat kernel 

p £ (x) = (2 7 re)- £i/2 exp(-|x| 2 /2e), 
in M. d . Then, we can consider the following family of random variables indexed 
by e > 

I £ (B H ,B H )^ f f p £ (B? - B?)dsdt, 
Jo Jo 

that we will call the approximated intersection local time of B H and B H . We 
are interested in the L 2 (fi) convergence of I £ (B H , B H ) as e tends to zero. 

For H = 1/2, the processes B H and B H are classical Brownian motions. 
The intersection local time of independent Brownian motions has been studied 
by several authors (see Wolpert [9] and Geman, Horowitz and Rosen [2J. The 
approach of these papers rely on the fact that the intersection local time of 
independent Brownian motions can be seen as the local time at zero of some 
Gaussian vector field. This approach easily allows to consider the intersection 
of k independent Wiener processes, k > 2. The applications of the intersection 
local time theory for Brownian motions range from the construction of relativis- 
tic quantum fields, see Wolpert [TU], to the construction of the self-intersection 
local time for the Brownian motion, see LeGall [4j. Further research has been 
done in order to study such problems for other types of stochastic processes, 
mainly Levy processes with a particular structure (strongly symmetric), see 
Marcus and Rosen [5J. 

In the general case, that is H ^ 1/2, only the self-intersection local time has 
been studied. Rosen studied in [TT] the planar case and a recent paper by Hu 
and Nualart [3] gives a complete picture for the multidimensional case. On the 
other hand, Nualart et al. [5] used a weighted version of the 3-dimensional self- 
intersection local time for the study of probabilistic models for vortex filaments 
based on the fractional Brownian motion . In recent years the fBm has become 
an object of intense study. A stochastic calculus with respect to this process 
has been developed by many authors, see Nualart [7] for an extensive account 
on this subject. Because of its interesting properties, such as short/long range 
dependence and selfsimilarity, the fBm it's being widely used in a variety of areas 
such finance, hydrology and telecommunications engineering, see |S]. Therefore, 
it seems interesting to study the intersection local time for this kind of processes. 

The aim of this paper is to prove the existence of the intersection local time 
of B H and B H , for an H ^ 1/2 and d > 2. We have obtained the following 
result. 
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Theorem 1 (i) If Hd < 2, then the family of random variables I £ (B H , B H ) 
converges in L 2 (f2). We will denote this limit by I(B H , B H ). 

(ii) If Hd > 2, then 

limE[I E (B H ,B H )} = +oo 

and 

\imV&T[I £ (B H ,B H )] = +oo. 

1 /2 

if {B; /A ,t>o} is a planar Brownian motion, then 




diverges almost sure, when e tends to zero. Varadhan, in [12j . proved that the 
rcnormalizcd self-intersection local time defined as lim £ ^o(I £ — exists in 

I? (Q). Condition (ii) implies that Varadhan renormalization does not converge 
in this case. 

For Hd > 2, according to the previous theorem, I £ (B H , B H ) doesn't con- 
verge in L 2 (17) and therefore I(B H , B H ), the intersection local time of B H and 
B H , doesn't exist. The proof of Theorem 1.1 rest on LemmalU which deals with 
the integral of a negative power of the determinant of some covariance matrix. 

The paper is organized as follows. In Section 2 we prove Theorem 1.1. In 
order to clarify the exposition, some technical lemmas needed in the proof are 
stated and proved in the Appendix. 

2 Intersection Local Time of B H and B , Case 

Hd < 2 

Let B H and B H two independent fractional Brownian motions on K d with the 
same Hurst parameter H £ (0, 1) . 

Using the following classical equality 

(2-Kj JR d 

from Fourier analysis, and the definition of I E (B H , B H ), we obtain 

I £ (B H ,B H ) =l[ T f T [ e^-^e-^dtdsdt. (1) 
(2tt) Jo Jo Jm d 

Therefore, 
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l[I £ {B H ,B H )] = -J-j f T f T f E[e^- B "- S ^}e-^^dsdt 
(2tt) Jo Jo Jk<* 

rp rp 

i / / / / I 2 tt | + 2/f> |£| 2 



(2n) a Jo Jo 

= 'r T r\e + s^ + t^)-^ d sdt, (2) 
(27Tj ' Jo Jo 

where we have used that (£, Bf - Bf ) - AT(0, |£| 2 {s 2H + i 2ff )), so 
E[e *(e,Bf-5f> ] = e -(.«+t«)i^ i 



and the fact that 



2H , .2H< [tlf / 27T 



d/2 



^e + s 2H + t 2H , 

According to the representation fl} for I e (B H , B H ), we have that 

E[I 2 (B H , B H )} — — ^— / / EljUtBf-SfWmBB-S?))] 
(2ir) J[o,T] 4 JR 2t! 

- i;i 2 +iii 2 ,,,,,,, 
x e 2 d^drjdsdtdudv. (3) 

Let introduce some notation that we will use throughout this paper, 

A = A(«,t) = s 2H + t 2H , 

p = p{u,v) = u 2H + v 2H , 

and 

/i = M (a, i, u, «) = \ (s 2H + t 2ff + u 2H + v 2H -\t- v\ 2H -\s- u\ 2H ) . 

Notice that A is the variance of B^ 1 ' 1 — B^ 1 ' 2 , p is the variance of Bff' 1 — B^> 2 
and p is the covariance between B^ 1 ' 1 — B^< 2 and B^' 1 — B^ :2 , where B H:1 and 
B H ' 2 are independent one-dimensional fractional Brownian motions with Hurst 
parameter H. 

Using that (Z,B? -B?) + ( V ,B? -B?) ~ N(0, A \tf + p\ v \ 2 + 2^,r,)) 
and © we can write for all e > 



E[I 2 (B H , B H 
1 



(27r) 2d J[o : T] 4 Jm.™ 



(2n) d J[o,T]' 



e -U(MeM 2 + ( P +e)M 2 +2^,v)} d ^ dr]dsdtdudv 

((A + e)(p + e) - t i 2 )- d/2 dsdtdudv. (4) 
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The last equality follows from the well known fact that 

e -i(w,Aw) dx 



(2n) d 



(det A) 
with 

'A + e p, 



1/2 ' 



A = IeL 



H p + e 



where Idd is the ci-dimensional identity matrix and £g> denotes the Kroncckcr 
product of matrices. We also have that 



det A = det Id d 



X + e p 
p, p + e 

A + e p 
p p + e 

((X + e)(p + e)-p 2 ) d . 



(det {Id d )) 2 ■ I det 



Proof of Theorem 1. Suppose first that Hd < 2. A slight extension of (f?|) 
yields 

E[I e (B H ,B H )I r ,(B H ,B H )] = [ ((\ + e){p + r 1 ) - fi 2 )- d/2 dsdtdudv. 

Consequently, a necessary and sufficient condition for the convergence in I? (fi) 
otI e (B H ,B H ) is that 

(Xp - p 2 y d l 2 dsdtdudv < +oo. 

[0 : T] 4 

Then the result follows from Lemma [4] 

Now suppose that Hd > 2, then from ([2]) and using monotone convergence 
theorem 

limni £ (B H ,B H )]= j T [ T (s 2H +t 2H y d/2 dsdt, 
J o Jo 

and this integral is divergent by Lemma [3] According to the expression (f2]) for 
E[I £ (B H ,B H )} and the expression © for E[I 2 (B H ,B H )} we obtain 

limVar[/ £ (B H ,i? H )] = lim S^E[I e {B H , B H ) 2 } - (e[I £ (B h , B H )} 

{Xp - p 2 y d/2 ~ [XpY d/2 dsdtdudv. 

[0,T] 4 

Set 

D £ := {(s, t, u, v) G R 4 + | s 2 + t 2 + u 2 + v 2 < e 2 }. (5) 
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We can find e > such that D e C [0, T] A . Making a change to spherical coor- 
dinates, as the integrand is always positive, we have 

/ (Xp - p 2 )- 1 ' 2 - (Xp)- 1 ' 2 dsdtdudv 

J[0,T] 4 

> f (Xp- p 2 )- d/2 - (M~ d/2 dsdtdudv = [ r 3 - 2Hd dr f * (6) d6, 

J D e JO J0 

where the integral in r is convergent if and only if Hd < 2, and the angular in- 
tegral is different from zero thanks to the positivity of the integrand. Therefore, 
ii Hd>2, then 

limVar[/ £ (B H ,B H )] = +oo. 

■ 

3 Appendix 

For clarity of exposition, we state and prove some technical lemmas in this 
appendix. 

Lemma 2 Let a > 0, and let 

l(a,x)= f e^y^dy (6) 
Jo 

be the lower incomplete gamma function. Then for all e < a and x > 0, 

7 (a, x) < K (a) x e , 
where K (a) = i V T (a) and T (a) = 7 (a, +00). 
Proof. If x > 1, 

7 {d, x) < T (a) x £ 7 
for all e > 0. On the other hand, if x < 1, 
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if £ < a. m 



(a,x)< / y ^dy= — = -x^ 
Jo a a 



Lemma 3 The following integral 



f f ( S 2H + t 2H y d/2 dsdt, 

Jo Jo 

is finite if and only if Hd < 2. 

Proof. It easily follows from a polar change of coordinates. 



G 



Lemma 4 Let 



{ {\p- n 2 )- d/2 dsdtdudv, 

J[0,T] 4 



then At < +00 if and only if Hd < 2. 

Proof. The necessary condition follows from a spherical change of coordinates. 
We can find e > such that D £ C [0, T] A , where D e is given in (0 . As the 
integrand in At is always positive we have 



where the integral in r is convergent if and only if Hd < 2, and the angular in- 
tegral is different from zero thanks to the positivity of the integrand. Therefore, 
if Hd > 2, then A T = +00. 

Suppose now that Hd < 2. By symmetry we have that 



T ee {(s, t,u,v) : < v < t,0 <t < T, < u < s,0 < s < T)}. 
Notice that 

A/9 - p 2 = detVar(Z) , 

where Z = (B?' 1 - Bf^ x ,B^' 1 - B"' 1 ). Due to the independence of B H and 
B H , we have that 





where 



Var (Z) = Var(Bf flf 1 ) + Var(Sf '\ B^ 1 ) 



and 



Xp - p 2 > det(Var(Bf B^ 1 )) + det(Var(flf ■ \ B* - 1 )), 

because the matrices Var(i?f r ' 1 . B^f' 1 ) and Var(_Bf r ' 1 , B^ 1 ) are strictly positive 
definite (see A8, (viii) in [5])- Then 




where 



ip (t, v) ee det(Var(Sf B* = i 2i V 



i(^+,--it-,r) 2 



Using Fubini's Theorem and 
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for all A, a > 0, we obtain 

(if (t, v) + <p (s, u))~ d / 2 dsdtdudv 
1 

1 

where 



r 

+00 

e -{v{t,-v)+<p(s,u))z z i-l dzdsdUudv 



T JO 

1 A 2 



zi- L A 2 (z)dz, (7) 







A(z) = ( [ e- v{t - v]z dvdt. 



'0 JO 

As A (z) < +oo, for all z S [0, 1], the integral © is convergent in a neighborhood 
of zero. Hence, we have to study the convergence of 

zi~ 1 A 2 (z) dz. 

Due to the homogeneity of order AH of tp (t, v) , if we make the change of coor- 
dinates t — z~4hx, v — z~4Hy, we obtain 

z*~ x A 2 (z)dz= J zJ-^ttIJ J e-^ x ^dydx\ dz. 

Now, using that {(x,y) : < x < Tz^n ,0 < y < x} C {(x, y) : x 2 + y 2 < 
2T 2 zzh} = S, and making a polar change of coordinates we have 

/ e-v^dydx < / e^^dydx = / re~ r " ] drd6 ', 

o Jo Js Jo Jo 

where <p (9) = f (cos 9, sin 9) . After the new change of variable x — r 4H (p(9) , 
the last integral is equal to 

,tt/4 x r 2 2H T 4H z V (6) , 

/ <p(0)~™ / —e- x xTm- 1 dxdd 
Jo Jo 4H 



tt/4 

V [ST™ 7 ((2H)" 1 , 2 2H T 4H z<p (6)) d9, 



1 

where 7 (a, cc) is given by ([S]) . Applying Lemma [51 
zi- x A 2 (z) dz 

1 /-+OC / pn/4 \ 2 

^ 7^72 / ( / P 7 ((2H)- 1 , 2 2ff T 4 ^ (0)) d0 ] 



16£P J 1 

16H 2 " \2H J J 1 



r)4Herp8He / i \ r+oo ( r -ir/i ^ 
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The integral in z is convergent provided e < 2 4 ^ d . It's an exercise of compu- 
tation of limits to prove that ip (9) ~ 8 2H as 9 [ and ip (6) - (tt/A - 9) 2H as 
9 | 7r/4, the main tool is to substitute the trigonometric functions by their first 
order approximations at the respective points. As a consequence, the integral 

ip(9) e ~™d6 
is always convergent. ■ 
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